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ABSTRACT 



We report new observations of CO (2-1) line emission toward five z~6 quasars 
using the Ka-band receiver system on the Expanded Very Large Array (EVLA). 
Strong detections were obtained in two of them, SDSS J092721. 82+200123. 7 and 
CFHQS J142952. 17+544717.6, and a marginal detection was obtained in another 
source, SDSS J084035.09+562419.9. Upper limits of the CO (2-1) line emission 
have been obtained for the other two objects. The CO (2-1) line detection in 
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J0927+2001, together with previous measurements of the CO (6-5) and (5-4) 
hues, reveals important constraints on the CO excitation in the central ~10 kpc 
region of the quasar host galaxy. The CO (2-1) line emission from J1429-I-5447 
is resolved into two distinct peaks separated by 1.2" (~ 6.9 kpc), indicating 
a possible gas-rich, major merging system, and the optical quasar position is 
consistent with the west peak. This result is in good agreement with the picture 
in which intense host galaxy star formation is coeval with rapid supermassive 
black hole accretion in the most distant universe. The two EVLA detections 
are ideal targets for further high- resolution imaging (e.g., with ALMA or EVLA 
observations) to study the gas distribution, dynamics, and SMBH-bulge mass 
relation in these earliest quasar-host galaxy systems. 

Subject headings: galaxies: evolution — galaxies: high-redshift — molecular data 
— galaxies: active — radio lines: galaxies 



Introduction 



Observations of high-redshift quasars probe the growth of supermassive black holes 
(SMBH) and their connection to galaxy formation at the earliest cosmic epochs. The dis- 
covery of strong sub millimeter /millimeter [(sub)mm] dust continuum in about 30% of the 
quasars known at z~6 provides the first evidenc e of active star forma tio n in young quasar 



host galaxies at the end of the reionization era (IBertoldi et al.ll2003al Jbl: iPetric et al.l 12003 



Priddev et al.ll2003l : iRobson et al.ll2004J : IWang et al. 



2007 



20081 ). The star formation rates 



estimated from the FIR luminosities (a few 10^^ to 10^^ Lq) are on the order of 10^ to 
lO^M0yr~^, which are comparable to the typical values found in so-called su bmillimeter 
galaxies at z = 2 ~ 3 fjScott et al.l l2002l : iGreve et al.l l2005l : iKovacs et al.l l2006l ) . The spa- 
tially resolved [C II] line emission from one of the most FIR luminous z~6 quasars, SDSS 
J114816. 644-525150. 3 (hereafter J1148-I-5251), further suggests a high star formation surface 
density of ~ 1000 Mr:, yr"^ kpc"^ o ver th e central 1.5 kpc region of the quasar host galaxy 
jMaiohno et al.ll2005[ Iwalter et al.ll2009l ). 



Mol ecular CO (6-5) line emission has been detected in ten of the FIR luminou s z~6 
quasars JBertoldi et al.ll2003bl : Iwalter et al.ll2003l : ICariUi eraDl2007l : IWang et al.l[201ol . 2011, 
in prep.), indicating the existence of highly-excited molecular gas in the quasar hosts. The 
CO (3-2), (6-5), and (7-6) transitions detected in the z=6.42 quasar J1148-I-5251 reveal 
a molecular gas component on scales of ~5 kpc in the host galaxy with CO excitation 
conditions simi lar to those found in local starburst galaxies and CO-detect ed quasars at 
lower redshifts JBertoldi et al.ll2003bl : IWalter et aP 120041 : iRiechers et aPboogh . 
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Emission in the low-order CO transitions (J < 2) from the z~6 quasar host galaxies 
is poorly con strained due to the limited sensitivi ty and frequency coverage of the previous 
instruments ( Wage et al. 20081 : Wang et al. 2010). T he new Ka band receivers on the Ex- 



panded Very Large Array (EVLA, iPerley et al. Il201ll ) open an imp ortant frequency window 



for studies of the co ld molecular gas in high-redshift galaxies (e.g., Ilvison et al.ll2010l . 12011 



Riechers et al.ll2010l ). In this paper, we report our EVLA observations of the CO (2-1) line 



emission in five z~6 quasars (iFan et al.l I2004J . l2006l: IWillott et al.l l2010al Jbr) Three of them 
are from the Sloan Digital Sky Survey (SDSS, Fan et al.l 12004 . 2006 ). with two objects, 
SDSS J084035.09+562419.9 and SDSS J092721.82+200123.7, previously detected in strong 
(> 3 mJy) 250 GHz dust continuum and molecular CO (6-5) and (5-4) line emission. Another 
object, SDSS J162331. 81-1-311200. 5, was detected in the [C II] line, but undetected in millime- 
ter dust continuum and high- J CO transitions (Bertoldi et al. 2011, i n prep.). The other tw o 
objects are from the Canada-France High-z Quasar Survey (CFHOS. IWillott et al.ll2010al Jbl) 
and do not have published CO observations yet. One of them, CFHQS 142952.17-1-544717.6, 
was detected in 250 GHz dust continuum (Omont et al. 2011, in prep.). We describe the 
observations in Section 2, present the results in Section 3, and discuss the CO excitation 
and host galaxy evolution properties of the detections in Section 4. A A-CDM cosmology 
with Hn = 71km s~ ^ Mpc~\ = 0.27 and Qa = 0.73 is adopted throughout this paper 
JSpergel et al.ll2007h . 



2. Observations 



The observations were carried out using the Ka-band receiver on the EVLA in 2010 in 
the D, DnC, and C configurations. The WIDAR correlator in Open Shared Risk Observing 
mode provided a maximum bandwidth of 128 MHz and a resolution of 2 MHz in each of 
the two basebands (A/C and B/D intermediate frequency [IF] bands). The A/C IFs could 
not be tuned below 32 GHz. The redshifts and observing frequencies of the CO (2-1) line 
of the five targets are estimated with previous detections of the CO (6-5), [C II], or quasar 
UV lines JCarilh et al.ll2007t IWang et al.lboioi Bertoldi et al 2011, in prep.). For the three 
sources with redshifts of z < 6.2 [corresponding to redshifted CO (2-1) line frequencies of 
t'obs > 32 GHz], we use the two 128 MHz IF pairs overlapped by 30 MHz and cover a total 
bandwidth of 226 MHz (i.e., ~2000 kms^^ in velocity and ~ 0.05 in redshift at z = 6)0. For 
the other two objects with z > 6.2, we centered the 128 MHz window of the B/D IF pairs 
on the line frequency and observed the continuum at >32 GHz with the other window. The 
observing time is 15 to 20 hours for each of the five targets (see Table 1). Flux calibrations 



^10 MHz overlap and a total bandwidth of 246 MHz for J1429+5447. 
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were performed using the standard VLA calibrators, 3C286 and 3C48, and we use 5-minute 
scan loops between targets and phase calibrators to calibrate the phase. The data were 
reduced with AIPS, and the spatial resolutions (FWHM) of the final images are typically 2" 
for data taken in the D configuration and 0.7" for the C configuration. 



Result 



CO (2-1) line emission has been detected in two of the five z~6 quasars, J0927+2001 
and J1429+5447, and marginally detected in J0840+5624. We present all the observing 
parameters and measurements in Table 1. The detailed results are listed below. 

J0927-I-2001 Toward this source strong dust continuu m at 850 GHz, 250 GHz, and 85 



GHz, and CO (6-5) and (5-4) line emission were detected (ICarilli et al.l 120071 : IWang et al. 



20101 ). We have detected the CO (2-1) line and the emission distribution (averaged over a 
velocity range of 880 kms~^) along with a spectrum is shown in Figure 2. The line peak 
emission centeroid is consistent with the optical quasar position and the peaks of the high-J 
CO lines. The line width (FWHM) and redshift fitted with a single Gaussian profile are 
590 ± 130kms~^ and 5.7716 ± 0.0012 which are in good agreement with the measurements 
from the high-ord er CO transitions (z = 5.7722 ± 0.0006 and FWHM = 600 ± 70kms-\ 
Carilli et al.l 120071 ). The line emission appears marginally resolved by the 2.19" x 1.96" 
synthesized beam with a peak surface brightness of 147 ±21 /iJybeam"^ and a total intensity 
of 230 ± 45 /iJy, with a source size of (2.7" ± 0.4") x (2.4" ± 0.3") determined from a fit with 
a two-dimensional Gaussian distribution (the deconvolved source size is about 1.7" x 1.4", or 
lOkpc X 8kpc). The corresponding line fluxes and luminosities (Table 1) are higher than the 
upper limits estimated from previous GET observations, but are still co nsistent given the 



large uncertaintie s and baseline feature contamination in the GET data (jWagg et al.l 12008 
Wang et ahlboioh . 



J1429-|-5 447 Toward this ob ject strong radio continuum em ission was detect ed in the 
FIRST survey (lEecker et al.lll995l ) and recent VLEI observations (IFrey et al.ll201ll ). making 
it the strongest radio source among the known z~6 quasars and the most distant radio-loud 
quasar. It has also been detected in dust continuum at 250 GHz with a flux density of 
~ 3 mJy (Omont et al. 2011, in prep.). We have detected both CO (2-1) line emission 
and continuum emission at the line frequency. The continuum source is unresolved by the 
0.71" X 0.67" synthesized beam and the flux density averaged over the line-free channels at 
32 GHz is 257 ± 15 /iJy. We subtract the continuum by performing linear fltting to the 
visibility data, using the UVLIN task in AIPS. The CO line emission is resolved into two 
peaks with a spatial separation of ~1.2" (6.9 kpc at the quasar redshift), and the optical and 
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radio quasar positions are consistent with the west peak (Figure 3). A Gaussian fit to the 
spectra yields a redshift of z = 6.1831 ± 0.0007 and a hue width of FWHM = 280 ± 70 km s"^ 
for the west source, and z = 6.1837 ± 0.0015 and FWHM = 400 ± 140 km s'^ for the 
east source. The hne fluxes estimated with the peak surface brightness on the velocity- 
averaged map averaging over a velocity range of ~ 450kms~^ are 0.065 ± 0.011 Jy km s~^ 
and 0.050 ±0.013 Jy kms~^ for the west and east components, respectively. However, a two- 
dimensional Gaussian distribution fitted to the east component suggest possible extension 
with a source size of (1.1" ± 0.2") x (0.7" ± 0.2"), which should be checked with deeper 
observations at higher spatial resolution. 

J0840-I-5624 This source was detected in (sub)mm dust continuum emission and CO 
(6-5) and (5-4) line emission; it has the broadest line width, FWHM = 860kms~^, among 



the CO-detected z~6 quasars ( IWang et al.ll2007l . |2010[ ). We observed the line at the redshift 
of z = 5.8441 ±0.0013 derived from the high-order CO detections and find no clear detection 
in a velocity-averaged map averaging over 1070 kms"^ made at the full resolution of 1.09" x 
0.76". At a lower resolution of 2.19" x 1.96", marginal signal (2.8cr) appears on the map 
(Figure 1), with a double-peaked mophology along the east-west direction. The optical 
quasar position is 0.8" away from the east peak. We plot the spectrum at the position of 
east peak in the right panel of Figure 1, and there is only very marginal signal (1 to 2a) over 
—500 to 500kni s~^, i.e., the typical velocity range of the CO (6-5) and (5-4) line emission 



( IWang et al.ll2010l ). The CO (2-1) line flux estimated with the surface brightness of the east 
peak is 0.062 ±0.022 Jykms"^ (Table 1). However, the signal is indeed marginal and deeper 
observations with a wider bandwidth are required to improve the measurement. 

J0210— 0456 This object is the highest redshift quasar kno wn to date with z = 6 .438 ± 
0.004 determined from the object's MgIIA2798A line emission Jwillott et al.lboiObh . We 
searched for CO (2-1) line emission in the 128 MHz window centered at the Mgll redshift 
but did not detect it. Here we assume a line width of 800kms~^, which is the typical full 
width at zero intensity (vfwzt) value found with samples of high-z CO-detected quasars 
( ICoppin et al.ll2008t IWang et al.ll2010[ ) to estimate the upper limit of the line intensity. The 
la rms noise level on the map averaged over this velocity range is arms = 16 /xJy beam~^, 
and the 3a upper limit of the line flux is estimated as 3crrmsVFwzi = 0.038 Jy km s~^. The 
corresponding 3a u pper limit of the line luminosity is L'co{2-i) < 1-28 x 10^° Kkms~^ pc^ 



(see equation (3) in lSolomon fc Vanden Boutll2005l ). However, we cannot rule out that the 
Mgll line emission is significantly offset from the quasar host galaxy redshift and CO (2-1) 
line falls outside the 128 MHz window. The continuum emission is also undetected with the 
other window centered at 32.1 GHz, and the channel-averaged map yields a 3a upper limit 
of < 54 yuJy. 
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J1623+3112 This object is detected in [CII] 158//m fine structure line emission by 



Berto ldi et al. (2011, in prep.), but undetected in 250 GHz dust continuum (IWang et al. 



20071 ). We searched for the CO (2-1) fine in the 128 MHz-bandwidth window centered at 



the [CII] redshift of z = 6.2605 ± 0.0005 and did not detect the fine. The rms on the map 
averaged over a velocity range of 800 km s~^ is 26 /iJy beam~^. This yields a 3a upper limit of 
< 0.062 Jykms"^ for the line fiux and < 2.0 x 10^° Kkms~^ pc^ for the line luminosity. The 
3a upper limit of the continuum emission at 35 GHz measured with A/C IFs is < 75 /xJy. 



4. Discussion 



We have observed molecular CO (2-1) line emission toward five quasars at z~6 using the 
EVLA, and detections/marginal detection have been obtained from the three objects that 
have strong FIR dust continuum emission. This is consistent with the picture of massive 
star formation fueling by huge amount of molecular gas in these young quasar hosts. The 
detection of [CII] in J1623+3112 is also likely to be a sign of star formation, but the current 
sensitivity of our EVLA observations cannot detect molecular CO from the host galaxy. 
J0927+2001 and J0840+5624, were previously detected strongly in the CO (6-5) and (5- 
4) transitions. CO (2-1) line emission has been detected and marginally resolved in the 
host galaxy of J0927+2001 over a scale of ~10 kpc. The molecular gas masses (Mgas) 
estimated from the CO (2-1) line peak surface brightness and the total intensity on the 
velocity- averaged map are listed in Table 1, assuming a CO luminosity-to-gas mass conversion 
factor of a = 0.8Mq (K km s~^ Pc^)~^ ap propriate for local ultraluminous infrared galaxies 



( ISolomon et al.lll997l : lDownes et al.lll998l ). These estimates are 1.7 and 2.5 time s higher than 



the v a lue of (1.8 ±0.3) x 10^° Mq estimated from the high-order CO transitions (jCarilli et al. 



20071 : IWang et al.l l2010l ). We plot the CO excitation ladder of this source in Figure 4, 
together with the results of Large Velocity Gradient (LVG) modeling of the highly-excited 
molecular gas components (gas densities of order lO^cm"^, kinetic temperatures of 50 to 
60 K, and peak at J > 6) found in other high-z FIR and CO lum inous quasars and nearby 
starburst galaxies (IRiechers et al.ll2006l . l2009l : iGiisten et al.ll2006l ). We normalize the models 
to the high-order CO transitions. The CO (2-1) line fiux measured with the peak surface 
brightness on the velocity-averaged map is consistent /marginally consistent with the values 
expected by these single-component models, while the total line fiux integrated over the 
line-emitting area falls above all the models. This may suggest the exsitence of additional 
low excitation gas i n the central ^ 10 kp c region as was found in the submillim eter galaxy 



AzTE C-3 at z=5.3 ( iRiechers et al.ll2010l ) and the nearby starburst galaxy M82 (IWeiB et al. 
20051 ). However, there are still large uncertainties in the measurements of all the three 



transitions, and observations of other CO transitions are necessary to address if there are 
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multiple CO excitation components in the quasar host galaxy. Our observations show no 
evidence of excess CO (2-1) hne emission and additional low excitation component in the 
host galaxy of J0840+5624. 

The C array imaging of the CO (2-1) line emission from J1429+5447 has resolved the 
molecular gas into two distinct peaks, with a spatial separation of ~6.9 kpc; the quasar 
position is consistent with the West peak. There is no clear velocity offset (26 ± 60kms~^) 
between the two components. These results suggest a gas-rich, major merging system with 
two distinct components that are comparable in CO luminosity and molecular gas mass. 
The west component of this system is in a radio-loud quasar phase. Similar quasar-starburst 
systems with multiple CO emission peaks were previously found i n the CO luminous quasars 
BRI 1202-0725 at z =4.7 JOmont et al.lll996l : ICarilh et al.ll2002h. B RI 1335-0417 at z=4.4 
f iRiechers et aPboOsI ) and J1148+5251 at z=6.42 dWalter et al1l2004h . These systems demon- 
strate the early phase of quasar-galaxy formation in which both AGN and starburst activities 
are tri ggered by major merger s and the molecular gas in the nuclear region is not fully coa- 
lesced ( iNarayanan et al.ll2008l ). We will expect further high- resolution observations with the 
EVLA in C or B array to constrain the gas surface density and dynamics, and with ALMA 
or the PdBI to resolve the dust continuum and distributed star formation in these young 
quasar host galaxies. 
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Table 1: Observations 



Name 


tobs 


Configuration 


Redshift 


FWHM 


lAv 


L'co(2-l) 


Mgas 




hour 








Jy kms~^ 


10l°Kkms-ipc2 


M© 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


SDSS J092721.82+200123.7 


20 


D, DnC 


5.7716±0.0012 


590±130 


0.129 ±0.018<= 
0.202±0.040'* 


3.70±0.52': 
5.80±1.15'' 


3.0±0.4= 
4.6±0.9'' 


CFHQS J142952.17+544717.6W'' 


15 


C 


6.1831±0.0007 


280±70 


0.065±0.011 


2.06±0.35 


1.6±0.3 


CFHQS J142952.17+544717.6E 


15 


c 


6.1837±0.0015 


400±140 


0.050±0.013 


1.59±0.41 


1.3±0.3 


SDSS J084035.09+562419.9 


20 


D, DnC, C 


5.8441'' 


860'' 


0.062 ± 0.022 


1.81±0.64 


1.4±0.5 


CFHQS J021013.19-045620.9 


15 


C 


6.438° 




<0.038 


<1.28 


<1.0 


SDSS J162331.81+311200.5 


20 


D 


6.26-'' 




< 0.062 


< 2.00 


<1.6 



"Redshift measured from the quasar Mgll line emission (Willott et al. 2010); ''Redshift and CO line width measured from the 



CO (6-5) and (5-4) Unes (Wang et al. 2010). '^CO (2-1) line flux and luminosity derived from the peak surface brightness on 
the velocity-averaged map; '*CO (2-1) line flux and luminosity derived from the total intensity integrated over the line-emitting 

area on the velocity-averaged map. '^The West component of J1429-I-5447 is consistent with the optical quasar position. 
■'^Redshift from the CII line detection (Bertoldi et al. 2011 in prep.) 
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Fig. 1. — Le/i(- Velocity-averaged map of the CO (2-1) line emission from J0840+5624 at a 
resolution of 2.19" x 1.96", averaged over —535 to 535 kms^^. The contour levels are (-2, 2, 
3)xl8 /iJybeam^^, and the la on the map is 21/iJy beam^^. The cross denotes the optical 
position of the quasar. Right-Spectrum at the position of the east peak, binned to a channel 
width of 20 MHz (178kms~^). The zero velocity corresponds to a redshift of z = 5.8441 
measured from the high-order CO emission (Wang et al. 2010). The error bar denotes the 
typical ±1(7 rms value. The dotted hue shows the zero intensity level and the sohd hue at 
the bottom denotes vpwzi of the high- J CO lines (Wang et al. 2010). 
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Fig. 2. — Le/t-Velocity-averaged map of the CO (2-1) line emission from J0927-I-2001 at a 
resolution of 2.19" x 1.96", averaged over ~ —440 to 440 kms""^. The contour levels are 
(-2, 2, 3, 4, 5)x21 /xJybeam"^. The cross denotes the optical quasar position. The source 
appears marginally resolved with a source size of 2.7" x 2.4" fit to a two-Dimensional Gaussian 
profile (deconvolved size of 1.7" x 1.4", or lOkpc x 8kpc). Right-Spectrum integrated over 
~2 beam area. The spectrum is binned to a channel width of 14 MHz (123 km s^^) and 
the zero velocity corresponds to the high-J CO redshift of z = 5.7722 (Carilli et al. 2007). 
The error bar represents the typical ±lcr rms in each channel. The dashed line denotes a 
Gaussian fit to the hue and the dotted hue shows the zero intensity level. The solid line at 
the bottom shows vpwzi of the high-J CO hues (CariUi et al. 2007). 
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Fig. 3.— Left-32 GHz continuum emission from J1429+5447 at a resolution of 0.71" x 0.67". 
The contour levels are (-2,2,4,8,16) x 15 //Jybeam"^. The cross denotes the optical position 
of the quasar. Middle- Velocity-averaged map of the CO (2-1) line emission, averaged from 
~ —100 to 340 kms^^ to get best S/N for both components. The contour levels are (-2, 
2, 3, 4)x27 /iJybeam"^. The source is resolved into two peaks with a separation of ~1.2" 
(6.9 kpc), and the quasar position is consistent with the west peak. Right-Spectra of the 
two components on the map. The quasar continuum has been subtracted from the West 
component. We adopt the surface brightnesses at the peak positions as the measurements 
of the line flux densities in each channel. The spectra are binned to a channel width of 14 
MHz (130kms~^) and the zero velocity corresponds to the redshift of z = 6.1831 measured 
with the west component. The dashed hues denote a Gaussian fit and the dotted fines show 
the zero intensity level. The spike at the 8th channel (i.e., ~ 500kms~^) of the upper panel 
is due to the poor sensitivity (lo" rms of ~ 160 fiJy beam^^) at the edge of the A/C IFs. The 
error bar represents the typical ±lcr rms in the rest channels. 



- 13 - 




Fig. 4.— CO excitation ladder of J0927+2001. The filled triangle at J=2 is the CO (2-1) 
line fiux measured with the peak surface brightness on the velocity-averaged map, while 
the open triangle is derived with the total intensity integrated over the line-emitting area 
of 2.7" X 2.4". The error bars show the la uncertainties of the line fiuxes. The solid and 
dashed lines show the LVG models of J1148+5251 with kinetic temperature, Tkin = 50 K, 
molecular hydrogen density, Pgas(H2) = lO'^'^ cm~^ from Riechers et al. (2009) and BR 1202- 
0725 with Tkin = 60 K, Pgas(H2) = 10"^'^ cm~^ from Riechers et al. (2006), respectively. The 
dotted line represents the LVG model of the central region of the nearby starburst galaxy 
NGC 253 with Tk^ = 60 K, Pgas(H2) = lO^-^ cm'^ from Glisten et al. (2006). These models 
are fitted to the CO (6-5) an (5-4) transitions. We have not made any correction for the 
continuum emission at the CO (2-1) line frequency as there is no significant detection. The 
continuum contimination to the CO (2-1) line fiux is estimated to be <5% based on the FIR 
dust continuum determined with previous (sub)mm detections (Wang et al. 2010). 



